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ON  THT-  MOTION  OF  ANCHOR  CHAINS 
3y  R.  Gran  Olsson,  Providence,  R.  I. 

1,  Introduction.  In  an  earlier  paper  the  problem  of  the  motion 
of  an  anchor  chain  was  treated  under  some  idealizing  conditions 
sue’1,  as,  for  instance,  neglecting  the  friction  of  the  anchor 
chain,  chain  holder  and  windlass  and  the  moment  of  inertia  of  the 
windlass  with  respect  to  the  rotating  axis  Cl],  The  anchor 
chain,  hov/ever,  also  moves  on  inclined  planes  with  different 
angles  and  fnc  ion  occurs  at  each  point  where  the  chain  is 
changing  its  direction.  It  may  be  shown  that  it  is  possible  to 
take  all  these  influences  into  account  without  making  use  of 
mathematical  tools  other  than  those  in  the  paper  cited  above.  In 
other  words:  it  is  possible  by  means  of  simple  substitutions  to 

consider  the  motion  of  anchor  chains  under  more  real  physical 
conditions  and  to  get  solutions  similar  to  those  obtained  earlier. 
The  results  may  be  interoreted  in  the  same  way  as  before. 

Bo  Tore  giving  the  mathematical  treatment  of  the  prob- 
lem let  us  consider  the  different  parts  of  the  anchor  mechanism 
and  the  conditions  vhic:'  this  mechanism  must  satisfy  in  order  to 
be  in  agreement  with  the  solution  given  in  this  paper.  The  linhs 
of  the  c ;ain  are  moving  from  the  chain  locker,  through  the  chain 
locker  pipe,  over  the  cha’-n  holder  of  the  windlass  and  through 
the  hawse  pipe  (Fig.  1)[2],  The  chains  should  have  a clear  lead 
to  the  hawse  pipe.  The  chain  holder  particularly  and  the  windlass 
generally  should  be  kept  well  lubricated  in  order  t,o  keep  the 
friction  moment  small.  The  position  of  the  chain  locker  is  very 
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important  if  the  windlass  gear  is  to  work  in  a satisfactory 
manner.  In(Fig.  2) is  shown  the  amount  of  grip  of  the  chain  Dy 
the  chain  holder:  to  insure  this  grip  it  is  essential  to  have 

the  tail  end  weight  of  the  chain  from  the  windlass  down  to  the 
chain  loclccr  well  below.  It  is  also  important  for  the  locker  to 
be  below  the  main  deck,  with  a chain  locker  pipe  (Fig.  1)  for 
controlling  the  lashing  cable  when  the  chain  is  dropped.  On 
merchant  marine  ships  anchor  chains  are  made  up  in  90  foot  lengths 
usually  coupled  together  by  a bolt  type  of  shackle  necessitating 
extreme-end  links  of  special  size  on  either  side  of  the  shackle. 
The  design  of  the  windlass  cable  holder  must  accommodate  the 
pitch  of  common  links  as  well  as  the  different  pitches  over  the 
shackle  connections  as  shown  in  (Fig.  2). 

The  transition  from  the  main  deck  to  the  hawse  pipe 
as  shown  in  (Fig.  3)  is  to  eliminate  the  wear  at  the  hav/se  pipe; 
to  reduce  the  frictional  force,  which  has  the  effect  of  decreasing 
the  effectivo  force  in  the  chain  between  the  ship  and  the  anchor, 
the  hawse  pipe  should  bo  given  such  a shape  as  to  provide  a 
smooth  and  continuous  race  for  the  chain  links. 


Notations . In  this  paper  the  same  notations  as  those  of 


Tim_oslien_ko  and  vqung  w ill  be  used  [3  ]•  some  more  mathematical 
notations  arc  in  agreement  with  those  used  in  textbooks  on  el- 
liptic integr;  Is  end  functions  [4  ], 
a = acceleration  (cm  sec  '), 
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C0  = constant  of  integration 
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e^,  (*2 » °3  e roots  of  the  equation  4-s^  - e 0. 

F » force  of  friction  (g), 

_2 

g = acceleration  of  gravity  (cm  sec  ), 

g2,  g^  = invariants  of  the  Weicrstrass  elliptic  function 

I = moment  of  inertia  of  windlass  (g  cm  sec  ), 

k = modulus  of  the  Jacobian  elliptic  function, 

= length  of  a section  o.f  the  chain  (cm), 

-1  2 

m = mass  of  the  system  in  motion  (g  cm  sec  ) 

2 

mred  c = rnass  °**  t*10  windlass  reduced  to  the  periphery  of 

-1  2 

the  windlass  (g  cm  sec  ), 
q =■  weight  per  unit  length  of  the  chain  (g  cm-*), 
r = average  value  of  the  windlass  radius  (cm), 

= tension  in  the  chain  for  an  arbitrary  section  (g), 
t time  (sec), 

u = parameter  of  the  elliptic  function  , 
v = velocity  (cm  sec"1), 

V - weight  of  the  anchor  (;,) 
o 

W = WQ  + gx  = weight  of  the  anchor  and  a length  x of  the  chain  (g), 
x = displacement,  length  of  the  chain  from  a fixed  point  (cm) 

2 

xQ  = WQ/g  = l°noth  of  the  chain,  corresponding  to  the  weight  W0(cm), 

= angle  of  an  arbitrary  inclined  plane 

Y = ratio  of  retarding  and  accelerating  forces 
A = discriminant  of  'weicrstrass  elliptic  function, 

5=  incriment  of  a quantity  (&x,  6q), 
e = angular  acceleration  of  the  windlass  (sec  ), 

= coefficient  of  friction, 

5 = (We  + gx  + g£  li  + gn»rcd)/W0, 
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9 = angle  of  rotation  of  the  windlass, 
w = the  real  half  period  of  the  elliptic  p-function 
p(u)  = the  elliptic  p-function  of  Veierstrass, 

C(u)  - the  elliptic  C-  function  of  Wciorstras3. 

3.  Equation  of  motion  and  its  integration.  By  using  d'Alemberts 
principle  and  the  principle  of  virtual  displacements  we  got  as 
equation  of  motion  of  the  system  by  an  increment  6x  of  displace- 
ment x and  6<p  of  angle  9 of  rotation  [5  ] . 

n 

(WQ  + gx  + g 2 ^ sin  a^)  6x  - F6x  - (I9  + Mp)  69 


- 7 4 [ (Vc  + qx  t qS^)  v ] 6x  = 0. 


(1) 


g <Tt 

As  will  be  shown  at  the  end  of  this  paper  it  is  not  necessary  to 
make  a nore  exact  consideration  of  the  friction  than  notating  it 
by  a constant  force  F. 

Because  of  the  fixed  connection  of  euchor  chair  and 
windlass  the  ’-incmatJc  condition  yields 


6x  = r&9 


(2) 


and  for  the  velocities  and  accelerations 

, « • • • • 

x = v 1*9  , x = a = r9, 

respectively. 

The  term  I9  69  may  be  written 

I x 6k  = m . 6x, 

rod  dt 


(2a) 


(3) 
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where  mrod  = I/r  is  the  moment  of  inortia  reduced  to  the  radius 
r of  the  windlass.  This  is  added  to  tho  last  terra  of  Eq.(l). 

In  order  to  obtain  a total  differential  the  expression 

Z?‘l±  (1  - sin  a1)  + gmred 

is  added  to  the  two  first  terms  of  Fq.(l),  giving  the  following 


equation 


w0  + qx  + qT'ti  + gmrod  - [F  + 4 (1  - sin  ai)  + gn^] 


i It  [(V'o  + + *Ui  + smTod)  v !• 


(la) 


! ultiplying  by 


(v.'o  + qx  + qlli  + gmrod)  v<H 


and  integrating,  we  obtain 


jg  (Wo  + q"-  + q^i  + b'mrod)3  - ^ [F  + ^ + q^i  (1  - sin  c^) 

+ Krared  ] (V/o  + ^ + °‘7\  + S»red  f * £g 

= {£:  (wo  + qx  + q^i  + gm  d)2,  (11 


where  the  integration  constant,  may  be  determined  from,  the  con- 
dition: v = 0 when  x - 0,  i.e, 

Co  = (’o  + qZti  + -nrcd)3  + C F + (1  - sin  a,) 

r 1 

+ Gmred  ] (Wo  + q7li  + Smred>2- 


J-L 
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It  is  convenient  to  introduce  the  following  abbreviations 


(W0  + qx  + + gm^)/^  * 


(4) 


<W0  + qS«t1  + gm^)/^  = K0l  (4a) 

[F  + + qS^d  - sin  a±)  + gni^  1/WQ  = Y.  (4b) 


The  velocity  may  now  be  written  as 


v2  = !l°  S 2 «3  _ S03)  + Y (5o3  - C2)Ta2. 
n 


With  the  notation 


Wc/g  = -o,  i.e.  1*  = d^, 


'o» 


*0 


(lc) 


we  obtain  after  some  simple  calculations  the  time 


t = Nj6  i f 


w 


c0  {Vu3  - C03)  + 6y  (C^T  t2")}17^  * 


(5)' 


when  it  is  assumed  that  £=  £0(~  = 0)  for  t = C. 

’’auction  (1)  is  now  integrated,  but  it  is  more  con- 
venient to  introduce  a new  independent  variable,  let  us  say  u, 
and  present  the  Mnoriatic  quantities  (displacement,  velocity, 
acceleration)  and  also  the  time  as  dependent  on  this  new  variable 
u. 


4.  Introduction  of  a nev;  Independent  variable.  From  Eq , (lb)  the 
time  t may  always  bo  expressed  by  the  integral 


-L  f < 


(6) 


V 

*‘0 
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where  x = xQ  correspond  to  t = 0.  In  the  second  integral  we 

have  introduced  C instead  of  x. 

Instead  of  expressing  the  time  as  a function  of  the  dis- 
placement let  us  introduce  a new  independent  variable 

„ = ['‘LtLalU  , (7) 


J i 


the  inverse  function  being 


C = p(u). 


From  Eq.(6)  we  have 


at  = a:  , 


and  from  Eq.(7) 

du  = &.  = 4£  = dt  . 

CC  C pf57 

Consequently  we  find,  by  integrating, 


t = 


tl 


u 


u. 


p(u)du. 


The  velocity  is  according  to  Eqs.(7a)  and  (0) 

C1  ^ = P'  (u)., 
du  p(u) 

The  acceleration  expressed  by  p(u)  is 


• • 

r _ <tt 
C " dl 


(8) 


(6a) 


(7a) 


(9) 


(10) 


•r  SL  (P/Xu))  = A.  p‘  (u)  du  _ p(u)p"(u)  - pV(u)2 
dt  p(u)  du  p(u)  dt  p(u)3 


As  C is  given  as  a function  F(C)  = F [p(u)]  we  get  the  following 
differential  equation  for  p(u)i 


UTS; 


8 


•iWL  '.''I  .- ■ 


_ , *~ry+  *= 


it: 


A 11 -91 


(ii) 


This  equation  may  serve  as  a cheek  at  the  end  of  the  calculation. 
In  our  case  v/e  have  to  substitute 


I 


l>(!^  - So5)  + 6y  (t?-  - to2)!’1' 


(7b) 


vhich  is  an  elliptic  integral,  and  the  inverse  function 


5=  P(u) 


(8a) 


is  an  elliptic  function  of  u. 
Py  means  of  the  substitution 


5 = (s  + ^),  d£  = ds, 

2 

we  transform  the  integral  (7b)  into 


(12) 


?T. 


^ 4 s^  - go  s - S3 


(7c) 


where  the  invariants  of  the  elliptic  integral 


’2  = 3/,  z3  - 0Sc^o  - | r>  * r3 


(13) 


and  the  inverse  function 


s = p(u) 


(8b) 


is  the  Weierstrass  p-function  if  we  choose  the  indicated  limits 
of  integration  in  (7b)  and  (7c). 


From  Fq.(5)  we  obtain  the  time 


|fr 


? 
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,«  - Y)  * t 5 ft ♦ i v2. 


By  means  of  Eq,(8b)  and  (12)  we  get 


a = dv  = 


4?  « — a — - t p(«)  - | y - l 

dt  p(u)  + J 2 6 (p(u)  * J)2 


(lib) 


(lie) 


There  is  an  asymptotic  value  of  acceleration  equal  to  g/3.  This 
constant  value  of  acceleration  corresponds  to  a straight  line  for 
velocity  and  a parabolic  curve  for  displacement. 

By  specializing  the  results  obtained  above  to  the 
case  Y = 0 we  get  the  same  expressions  as  in  the  earlier  paper. 
Therefore  we  don't  need  consider  this  special  case  now. 

5.  jAAsA4aA44 °A  A4  Ah.e.  d.l 3.c„r iroiASAt  jc£_Jthe_  We i e rsJ(.rAsA  P- 
f unction.  From  a numerical  calculating  point  of  view  it  is  im- 
portant to  consider  the  discrimant  of  the  1 cierstrass  p-function, 
i.e.  [ 6 ] 

& = s23  - 27s32,  (14) 

where  g2,  are  the  invariants,  introduced  in  Bq.(7c).  Insert- 
ing g2  and  g^  from  Eq.(13)  in  (14)  we  obtain  after  some  elemen- 
tary calculations 

A = - 216  502(*;0  - j Y)U3  + 2 - 3 Y ^2),  (14a) 

where  ^ and  Y are  given  by  (4a)  and  (4b),  From  Eq.(la),  however, 
we  see  that  motion  is  nossible  only  if  K0  > Y and  from  (4a)  and 
(4b)  that  > 0,  Yj>0,  where  Y = 0 corresponds  to  the  case  in- 
vestigated earlier.  Writing  the  discriminant 
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a = ^32  5o6c|  .£  - 1)<£  - 1)2(^L  + 1).  (14b) 

The  value  of  the  discrimant  for  Y = 0 is 

4o  " - ^2  c06.  (14c) 

In  order  to  obtain  a convenient  representation  of  the  discriminant 
A we  write 

A/^o  = " (2  f*  “ + 1)*  (l1*!) 

This  expression  considered  as  a function  of  X/K0  has  a maximum 
for  - ( \fl?  - 1)A  = 0,?8l  and  is  equal  to  zero  for 

Y/^o  “ 2^3  anci  Y^0  ~ The  curve  A/Aq  is  shown  in  Fig.  5 . 

Ix  we  write  the  polynomial  under  the  square  root  of 

Eq. (7c) 


- Z23  ~ S3  = **(s  - e1)(s  - e2)(s  - e,)*  (15) 

el»  e2 » e3  beinS  the  roots  of  the  equation 

kz  - g2s  " S3  = 0 (15a) 

the  negative  value  of  the  discrimantA  means  that  two  of  the  roots 
let  us  say,  el  and  are  complex,  the  third  root  e?  being  real. 
Introducing  the  quantities  H and  k2  defined  by 


(ID 


(15a) 


V‘  = c2  - e1)<ea  - c3)  . 2e22  ♦ 

k2  = 1 - 3e2 

2 i?sr» 


(16) 
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where  k Is  the  modulus  of  the  Jacobian  elliptic  functions,  Eqs. 
(16)  establish  the  connections  between  the  quantities  g^,  ep,  H 
and  k2,  where  e 3Y^  and  e0  may  be  found  as  the  real  root  of 
Eq.(l5a).  To  get  a survey  over  the  numerical  calculations  it  is, 
however,  more  convenient  to  consider  the  values  H and'k2  as  given 
end  from  these  calculate  g^  and  e2,  getting 


e2  = | H(1  - 2k2), 

g3  = (2H)3(1  - 2k2)  [9  - 8(1  - 2k2)2]  . 


(16a) 


The  connection  between  the  invariant  g2  and  g^,  on  the 
other  hand,  are  given  by 

g2  « lte22  - fa  = 3y2.  (17) 

e2 

Introducing  e2  and  g^  from  ^qs.(l6a),  we  obtain 

S2  • ‘*H2[  | (1  - 2k2)  - 1]  . (17a) 

Because  of  g2  = 3Y2,  we  must  have  g2  > 0 or 
| (1  - 2k2)2  - 1>  0, 

which  neans  k ^ 0.06699,  k < 0.259,  the  sign  equal  corresponding 
to  the  equianharmonic  case  (g2  = 0)  considered  earlier.  Because 
of  g2  >0  the  value  of  k2  is  3mall  and  therefore  as  a good  ap- 
proximation we  may  assume  1c2  = 0,  giving 


t 

m 
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52  « | H = ^g2  « Y, 

a - (S  U*»-  - y3 
u3  ^5  r* 


On  the  other  hand  we  have  from  Eq.(13) 


(16b) 


g3  -«0  «o  “2Y)  + Y3 


(13a) 


and  therefore  the  following  condition  between  £Q  and  Y holds  for 
the  case  k“  = 0: 


= 2/3, 


all  quantities  being  expressed  by  Y. 


As  Y/£  = 2/3  corresponds  to  the  special  case  A = 0, 


where 


4s3  - g2s  - =•  4S3  - 3Y2s  - Y3  = 0 

with  all  roots  real  e1  = y,  e2  = = - Y/2,  we  have  the  follow, 

ing  equation  for  the  Weierstrass  p-function 


p(u)  = e^  + — eJ.  Z ,e2 

sn2(u  fa  - e^) 


(18) 


the  real  half  period  being 


Ut  = 


2L 


{ “ r3  v*  ^Y  \J  bej_ 


(19) 


6 . Representation  of  the  ; -inenatlc  quantities  in  the  case 

VL  = 2/3. 
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From  Eq.(l8)  we  get,  Inserting  for  and  e^ 

Y £ Y 

p(u)  ♦ - ' • 

- sir*  (u  w ^ T ) 


The  value  of  u corresponding  to  the  argument  V2  is  u0  = 


and  therefore  wo  may  write 


J\ T 


p(u)  + * = 


2 sl"2<f  & 


Introducing  this  in  Eq.(Cc)  we  get,  remembering  from  Eq.(4a) 
that 


Lli  + q mred  " — "*o 


- ^ «0  - 1)  - * - 1). 


(8d) 


x ^ x0  i Y [ 1 

0 2 o r* — n 

sin2(u  \j  i Y ) 


- 1 1 


From  this  ve  see  that  the  boundary  condition  x = 0 for  t = 0 


(u  = uQ)  is  satisfied,  Equation  (Od)  may  be  written 


x = 


Y x0  cot2  (u  >||  Y ). 


(8a) 


The  time  may  bo  written  as 

t . 4°,^  f t p(u,  + x , 

s i 2 


du  » C^2)1/2  3 Y 
5 2 \] 


ft 


du 


sin 


( \J|  ? u) 


(5d) 


which,  after  integration,  yields 


L 
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t = - ctis)1/2  (|  r )1/2  cot  <4§7u). 


From  this  we  obtain 


^ Y cot2(  vl  5 Y u)  c 


(5e) 


(5f) 


which  inserted  in  Eq.  (Oe)  yields  the  simple  result 


x = | t . 


(8f) 


From  Eq.(lOa)  for  the  velocity  we  get 


x = - 2 (2  r)1/2  cot  JJ7 u = 


2(-6^)1/2  (iX-2)“1/2t 


(10b) 


or,  in  agreement  with  (Of), 


x = fi  t. 
3 


( 10c ) 


In  the  same  manner  v/e  finally  obtain  for  the  acceleration 

cos2 ( \||  y' u) 


a = ^ 5 sin2(  v(|  Y u) 


L sin 


(JJ7  u 


- Y -Y 


a - 2 2 sin2( 
3 Y 


Y u)  [ ^ y - Y ] 


sin2(  ^ Y* u) 


sin2(\|  3 Y^ u 


I’ 


) 


V/e  therefore  have  for  Z,Q  = | y the  very  simple  case  of  constant 
acceleration,  velocity  linearly  dependent  on  time  and  displace- 
ment depending  on  the  square  of  time.  If  we  are  starting  with  an 
acceleration  equal  to  g/3,  the  motion  will  continue  with  this  con- 
stant value  of  acceleration.  In  Fig.  6 is  3hown  the  displacement 
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as  a function  of  time  for  the  cases  Y/<0  B 0,  Y/^0  • 1/3  and 
Y/50  = 2/3  (constant  acceleration) $ Fig.  ? and  8 present  velocity 
resp.  acceleration  for  the  same  values  of  Y/£0# 

7.  Simple  proof  that  constant  acceleration  Is  .equal  _to_jg/^.  The 
equation  of  notion  may  be  written 

4fv  + m ft.P,  (18) 

where  dm  = qdx/gt  dra/d t = qv/g, 

dm  v = 3 v2.  (19) 

dt  g 


At  the  beginning  of  motion  we  have  v ■ 0 and 


dv  c dv  _ o 

at  an  ' 

(t=0) 


5t=o  * m°  a?  “ “oa- 
(t=0) 


(18a) 


Substracting  this  last  equation  from  Eq.(lS)  we  get 


a v + (m  - mQ)  a = P - P,  . = qx  (l8b) 

g ° (t«sO) 

because  the  difference  of  forces  at  an  arbitrary  time  t and  t = 0 

is  qx.  Furthermore  is  m - mQ  = 2 x.  By  constant  acceleration  a 

? 3 

the  square  of  velocity  is  v = 2ax,  which  value  inserted  in  Eq. 
(l8b)  yields 


- 2ax  + SI  ax  « qx  or  a 
Z 8 


= const. 


This  is  the  only  value  of  constant  acceleration  which  is  valid 
for  all  values  of  x,  that  means  during  the  whole  motion  from 
beginning  to  stop.  From  Eq.(l8b)  we  see  that  a constant  acceler- 
ation i3  possible  by  motion  with  variable  mass  if  the  mass  in 
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motion  and  acting  forces  are  increasing  proportional  to  the  dis- 
placement. If  the  notion  takes  place  under  the  influence  of 
gravity  the  constant  value  of  acceleration  is  just  equal  to  g/3. 

8.  Summary.  The  motion  of  an  anchor  chain  under  the  influence 
of  gravity  and  friction  forces  is  considered.  The  equation  of 
motion  involves  also  the  influence  of  the  moment  of  inertia  of 
the  windlass  and  the  frictional  moment  acting  on  the  axis  of  the 
windlass.  The  solution  of  the  equation  can  be  expressed  by  the 
elliptic  functions  of  Weierstrass,  the  time  by  the  C-  function 

of  a parameter  u,  the  kinematic  quantities  (displacement,  velocity, 

acceleration)  by  the  p-function  and  its  derivative  of  the  same 

parameter.  In  a special  case  these  functions  can  be  expressed  by 

trigonometric  functions,  corresponding  to  a constant  acceleration 

2 

a = g/3,  the  velocity  v » gt/3  and  the  displacement  s = gt  /6, 

The  absence  of  frictional  forces  corresponds  to  the  equianharmonic 
case  of  the  Weierstrass  function,  tabulated  by  A.  G.  Grcenhill 
[7  3.  All  real  cases  of  motion  of  the  anchor  chain  are  included 
between  this  case  and  the  case  of  constant  acceleration  a = g/3 
and  can  be  found  by  linear  interpolation  with  respect  to  the 
parameter  Y,  expressing  the  ratio  of  the  friction  forces  and  the 
accelerating  forces. 

9.  Appendix.  As  Mentioned  in  section  3 it  is  possible  to  give 
a more  complete  expression  of  the  frictional  force.  This  will 
be  done  in  the  following.  Let  us  consider  the  forces  acting  on 
an  inclined  plane  (n  - k)  with  the  angle  of  the  slope  equal  to 
a(n  k)’  tV>:C  l°ngth^(n  ^ and  the  friction  coefficient  p , which 


* . 
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may  be  assumed  equal  for  all  inclined  planes  (Fig, 9).  For  the 
motion  of  this  part  of  the  chain  we  have  the  followinc  equation 

sn-k  - sn-k  + e^-n-k  (sln  Vk  * •*  003  °n-k)  " § *h-k  If-  (20) 


At  the  corner  of  two  adjacent  slopos  there  exists  the  following 
relation  between  the  forces  in  the  chain  on  both  sides  of  the 
corner 


S'  . = S exp  [n(a 

n-k+i  n-k  n-k+1 


- an-k>  ] • 


(21) 


The  prime  index  indicates  the  force  in  the  chain  just 
below  the  corner,  the  same  notation  without  the  prime  index  in- 
dicates the  force  on  the  same  slope  just  above  the  corner.  If 
the  first  inclined  plane  where  the  chain  is  moving  after  leaving 
the  windlass  has  the  number  (n  - m),  the  force  in  the  chain  is 

SU  ■ v)  + £ * lu°  + *1  + ' 

> (22) 

- q(^Q  sin  a0  + ^1  sin  a1  + ^2  sin  a2^* 

^o*  avC:  the  c^ain  the  three  first  sec- 

tions, according  to  Fig.  4.  On  the  other  hand  the  force  in  the 
chain  just  below  the  last  corner  corresponding  to  the  vertical 
plane  (n  + 1),  is  given  by 

sn+l  = (Wo  * 9X)  - | ^ Kwo  + <ix)v  ] . (23) 


4 * 


From  the  last  four  equations  we  can  establish  an  equation  of 
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motion  similar  to  Sq.(la),  section  3.  It  is  convenient  to  in- 
troduce 


L 


™ r (| 


a 

n-m 


>3 


L' 


k-m 

+ T l . exp 

!c=0  n K 


k=p 

T 

ktrO 


'(n-k) 


(sin  a 


n-k 


- H 


t*1'?  ■ a(n-k)  3) 
cos  an-k)  exp  Cn(g 


U24) 


If  wc  expand  the  equation  (23)  on  both  sides  by  (guty^  ♦ qL),  in 
order  to  have  n total  differential  we  obtain  the  equation  of 
motion 

(WQ  + qx  ♦ gmred  + qL)  - i ^ ((w0  + qx  + gmred  + qL)  v] 

= emred  + ’L  ~ {r  eXp  ftl<i  ‘ “n-m5  3 * <?L’  ]'•  C25) 

The  two  last  terms  contain  the  influence  of  the  friction  forces 
on  the  system.  If  we  now  multiply  the  equation  by 

(Wq  + qx  + gmred  * qL)  vdt  = (WQ  + qx  + Z*Ted  + qL)  dx 

and  integrate,  we  get  the  same  equation  as  Eq.(lb),  section  3» 
apart  a slight  change  of  some  of  the  constants.  We  thus  have 
given  a proof  that  it  is  allowed  to  introduce  the  friction  as  a 
constant,  if  we  neglect  the  influence  of  centrifugal  forces  at 
the  corners  during  motion. 
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